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ABSTRACT 


A  method  is  set  forth  for  the  computation  of  test-section,  after¬ 
shock,  and  stagnation  conditions  in  hotshot  tunnels  of  Arnold  Engineer¬ 
ing  Development  Center,  using  nitrogen  at  stagnation  temperatures 
from  1500  to  5000°K.  The  basic  input  data  are  the  stagnation  pressure 
in  the  reservoir,  the  pitot  pressure  behind  a  normal  shock,  and  the 
heat -transfer  rate  at  the  stagnation  point.  The  value  of  the  overall 
enthalpy  is  adjusted  by  numerical  iteration  until  the  measured  and  com¬ 
puted  heat-transfer  rates  are  numerically  equal  within  a  prescribed 
tolerance. 
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NOMENCLATURE 

Effective  area  ratio 
Speed  of  sound,  ft/sec 

Parameters  in  empirical  equations  in  Appendix  I,  with 
numerical  subscripts 

Throat  diameter,  in. 

Enthalpy,  atm/amagat  units 

Enthalpy,  ft 2 /sec 2 

Tunnel  constant 

Dimensionless  enthalpy  function,  hp/p 

Logarithm  to  base  10 

Mach  number 

Symbol  for  the  ratio  k/yp; 

Dimensionless  parameter,  cf.  Eq,  (9) 

Pressure,  atm 

Heat-transfer  rate,  Btu/ft^sec 

Gas  constant,  equal  to  3.  661  (10)~2  when  pressure  is  given 
in  atm,  density  in  amagats,  and  temperature  in  degrees 
Kelvin 

Radius  of  model,  in. 

Reynolds  number  per  foot 

Density  ratio  across  normal  shock,  P^/Pg 

Entropy,  dimensionless 

Temperature,  °K 

Longitudinal  flow  velocity  in  units  such  that  u2  is  given  in 
atm/amagat  units 

Longitudinal  flow  velocity,  ft /sec 

Average  longitudinal  signal  flow  velocity,  ft/sec 

Symbol  for  the  quantity  (1/4)M2^ 

Ratio  of  specific  heats 

Isentropic  expansion  coefficient 
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p  Viscosity,  lb/ft-sec 

p  Density  in  amagat  units,  a  relative  scale  referred  to 

density  at  0°C  and  one  atm  pressure.  This  scale  is  the 
same  as  the  ratio  p/pQ  in  Refs.  6  and  7,  using 
pQ  =  0.  00242  Ib-sec2/ff4  for  nitrogen. 

SUBSCRIPTS 

o  Relates  to  stagnation  conditions.  On  an  unprimed  quan¬ 

tity  (e.  g, ,  pQ)  it  denotes  reservoir  conditions.  On  a 
primed  quantity  (e.  g.  ,  pQ')  it  denotes  stagnation  conditions 
behind  a  normal  shock. 

1  Condition  just  before  a  normal  shock 

2  Condition  just  behind  a  normal  shock 

Numerical  subscripts  are  also  used  to  identify  param¬ 
eters  bj  and  c^  in  empirical  equations  in  Appendix  I. 


SUPERSCRIPT 


Throat  condition 
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1.0  INTRODUCTION 


Based  in  part  on  previous  work  (Ref.  1),  this  report  presents  an 
improved  method  of  determining  the  test-section,  after-shock,  and 
stagnation  conditions  in  hotshot  tunnels  at  the  Arnold  Engineering 
Development  Center  (Ref.  2),  using  nitrogen  as  working  gas  at  stagna¬ 
tion  temperatures  from  1500  to  5000°K.  Besides  stipulating  thermo¬ 
dynamic  equilibrium  at  all  times,  as  is  also  being  done  here,  the  earlier 
paper  assumed  that  no  temperature  or  density  gradients  exist  in  the  arc 
chamber  and  that  the  decay  of  reservoir  conditions  as  a  function  of  time 
can  be  calculated  from  the  rate  of  mass  flow  in  the  throat. 

These  latter  assumptions  became  suspect  in  a  roundabout  way 
(Ref.  3).  In  run  after  run,  the  results  exhibited  large  and  erratic  dis¬ 
crepancies  between  the  measured  and  computed  rates  of  heat -transfer 
at  the  stagnation  point,  the  calculated  values  being  derived  from  Fay- 
Riddell  theory  (Ref.  4).  At  first  these  discrepancies  were  attributed 
solely  to  contamination  of  the  working  gas  caused  by  erosion  of  metallic 
surfaces  in  the  arc  chamber.  As  is  described  in  Ref.  3,  successive 
new  designs  of  the  arc  chamber  diminished  the  contamination  and  also 
reduced  the  differences  between  the  measured  and  calculated  values  of 
the  rates  of  heat  transfer  at  the  stagnation  point.  But  as  further  new 
designs  of  the  arc  chamber  reduced  the  contamination  substantially, 
residual  discrepancies  remained  between  the  measured  and  computed 
rates  of  heat  transfer,  which  could  no  longer  be  ascribed  to  contamina¬ 
tion. 


Suspicion  now  fell  on  the  supposed  uniformity  of  conditions  in  the 
arc  chamber,  and  at  this  juncture  it  appeared  advisable  to  measure  one 
or  another  parameter  in  the  free  stream.  This  was  successfully  done 
(Ref.  3)  by  means  of  timed  photographic  records  of  small  blast  waves  in 
the. free  stream,  which  were  shown  to  be  moving  with  the  local  flow 
velocity.  However,  this  flow  velocity  in  the  free  stream  could  also  be 
calculated  in  two  ways,  one  in  terms  of  the  supposedly  uniform  reservoir 
conditions  and  the  other  in  terms  of  the  measured  rate  of  heat  transfer 
at  the  stagnation  point.  The  values  derived  from  the  heat -transfer  rate 
consistently  agreed  with  the  measured  velocities  to  within  5  percent, 
whereas  the  values  based  on  reservoir  conditions  sometimes  differed 
from  the  measurements  by  as  much  as  20  percent. 

This  led  to  the  conclusion  that  the  measured  heat -transfer  rate  at 
the  stagnation  point  could  be  adopted  as  a  basic  input  datum,  which 
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would  determine  the  system  if,  in  addition  to  the  measured  pitot  pres¬ 
sure  downstream,  the  measurement  of  one  reservoir  condition  could  be 
accepted.  By  virtue  of  the  high  speed  of  sound  in  the  reservoir,  it  is 
reasonable  to  assume  that  the  pressure  gradients  in  the  reservoir  dis¬ 
appear  in  a  few  milliseconds.  The  flow  in  the  tunnel  is  not  fully  estab¬ 
lished  until  about  30  milliseconds  after  the  arc  is  fired,  so  the  pressure 
measurements  are  deemed  to  be  trustworthy.  However,  heat-transfer 
theory  suggests  that  the  temperature  gradients  might  endure  for  one 
second  or  longer,  so  the  density  gradients  would  be  equally  slow  in  dis¬ 
appearing. 

The  indicated  procedure  is  then  as  follows:  Fay-Riddell  theory 
(Ref.  4)  determines  a  value  of  the  stagnation  enthalpy  based  on  the 
measured  heat-transfer  rate  and  pressure  at  the  stagnation  point. 
Corresponding  values  of  density,  temperature,  and  entropy  in  the 
reservoir  are  then  given  by  empirical  equations  written  in  terms  of 
pressure  and  enthalpy  as  independent  variables.  Thereafter,  a  com¬ 
putational  program  not  unlike  the  one  of  Ref.  1  determines  the  system, 
including  a  calculated  value  of  the  heat-transfer  rate  at  the  stagnation 
point.  If  this  computed  value  of  the  heat -transfer  rate  is  unacceptably 
different  from  the  measured  value,  the  computer  is  instructed  to  vary 
the  stagnation  enthalpy  until  the  measured  and  computed  values  of  the 
heat -transfer  rate  are  in  agreement  to  within,  say,  1  percent.  Inas¬ 
much  as  the  ultimate  value  of  the  overall  enthalpy  is  found  by  iteration, 
the  value  of  enthalpy  used  to  start  the  iteration  need  be  only  approxi¬ 
mately  correct.  This  procedure  is  independent  of  any  temperature  or 
density  gradients  in  the  reservoir.  It  also  renders  it  unnecessary  to 
compute  rates  of  decay  of  reservoir  conditions  because  instantaneous 
values  of  these  conditions  are  derived  from  input  data.  However,  in 
selecting  values  of  pressure  in  the  reservoir,  allowance  must  be  made 
for  the  time  of  passage  of  a  pressure  signal  from  the  throat  to  the  stag¬ 
nation  point  downstream  (see  Ref.  1). 

Data  for  the  thermodynamic  properties  of  nitrogen  are  taken  from 
the  tables  of  Smith  (Ref.  5),  Hilsenrath  and  Klein  (Ref.  6),  and  Little 
and  Neel  (Ref.  7).  The  tables  of  Smith  are  based  on  the  virial  coeffi¬ 
cients  given  by  Amdur  and  Mason  (Ref.  8),  and  they  do  not  account  for 
the  effects  of  dissociation  and  ionization.  The  referenced  tables  of 
Hilsenrath  and  Klein  account  for  intermolecular  forces,  as  well  as  dis¬ 
sociation  and  ionization,  but  their  range  of  densities  extends  upward 
only  to  about  200  amagats. 

The  computational  program  of  this  report  is  intended  for  the  follow¬ 
ing  ranges  of  conditions:  The  temperature  in  the  reservoir  after  firing 
the  arc  is  assumed  to  be  between  1500  and  5000°K,  and  the  corresponding 
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range  of  pressures  extends  from  10  to' 2500  atm.  Upstream  of  the 
shock,  the  temperature  and  density  are  assumed  to  be  between  the 
liquefaction  point  and  400°K  and  between  10~5  and  10"  1  amagat,  respec¬ 
tively.  Behind  the.  shock  the  pressure  is  assumed  to  be  between  10“^ 
and  one  atm. 

A  detailed  outline  of  the  computational  program  is  given  in  Appen¬ 
dix  I,  In  order  not  unduly  to  encumber  the  ensuing  text,  some  of  the 
equations  in  the  program  are  given  only  in  the  appendix,  with  appropriate 
references  in  the  text.  Furthermore,  lest  certain  equations  in  the  text 
be  littered  with  conversion  factors,  enthalpy  (h)  and  the  square  of  the  flow 
velocity  (u2)  are  often  expressed  in  atm/ amagat  units._  However,  in  the 
working  equations  in  Appendix  I,  enthalpy  appears  as  h  and  the  square  of 
the  flow  velocity  as  u  ,  both  being  expressed  in  ft2/sec2.  Some  of  the 
empirical  equations  in  Appendix  I  were  written  by  means  of  the  method  of 
surface-fits  developed  by  Lewis  and  Burgess  (Ref.  9).  Appendix  II  com¬ 
prises  computer  printouts  from  a  number  of  typical  runs  and  for  a  variety 
of  input  conditions. 


2.0  STAGNATION  POINT  HEAT-TRANSFER  RATE  AND 
RESERVOIR  CONDITIONS 


In  the  theory  of  Fay  and  Riddell  (Ref.  4),  the  boundary-layer  equa¬ 
tions  are  developed  for  the  case  of  high-speed  flight,  in  which  the  fluid 
in  the  external  flow  is  partially  dissociated.  The  analysis  includes  the 
effects  of  diffusion  and  recombination,  and  it  finally  reduces  to  a  set  of 
non-linear  equations  which  are  amenable  to  numerical  solution.  If  the 
effects  of  dissociation  are  neglected,  and  if  the  Prandtl  number  is 
assumed  equal  to  0.  71,  the  rate  of  heat  transfer  at  the  stagnation  point 
is  given  by  the  relation 


q  VRa  =  constant  x  (pwtUv  )0,1  (poX)  (h0  -  hw) 

,  ,  .0.25  ,  ,v~0.25 

x  (,p0  -  p,)  Kpo) 


(1) 


in  which  q  is  the  heat -transfer  rate  and  Ra  is  the  radius  of  curvature  of 
the  model  surrounding  the  heat-transfer  gage.  The  temperature  at  the 
wall  is  assumed  to  be  300°K,  and  the  pitot  pressure  is  also  the  wall 
pressure.  Under  these  conditions,  the  density  at  the  wall  is  proportional 
to  the  pitot  pressure,  whereas  the  viscosity  and  enthalpy  at  the  wall  are 
constants.  Equation  (1)  then  becomes 


q  V  Ra  =  constant  x  (uo')0,<  (p o)°'*5  (h0  -  hw) 
X  (po)0'1  (po  -  Pl)°'a5 


(2) 
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In  this  equation,  if  q  is  expressed  in  Btu/ft^sec,  Ra  in  inches,  p  in 
amagats,  M  in  lb/ft-sec,  p  in  atm,  and  h  in  ftS/sec^,  the  constant  is 
numerically  equal  to  4.  2519(10)"4  and  the  enthalpy  at  the  wall  is  equal 
to  3.  3469(10)6.  In  the  case  of  nitrogen,  the  conversion  from  h  (given 
in  atm/amagat  units)  toll  (ft^/sec^)  follows  from  the  relation 

h  =  8.722  (  10)5  h 

In  the  range  of  stagnation  temperatures  from  3000  to  4000°K, 
numerical  evaluations  of  Eq.  (2)  in  more  than  50  runs  disclose  the 
empirical  fact  that  the  stagnation  enthalpy  is  approximately  a  linear 

function  of  the  quantity  (q  V  Ra/pQ'  ).  Accordingly,  a  provisional  value 
of  the  stagnation  enthalpy  is  given  by 

ho  =  1.459  (10')*  (q  V  Ra/po)  +  0.7750  (10)7  (3) 

By  using  this  value  of  hG  and  the  measured  value  of  pQ  in  the  reservoir 
as  independent  variables,  the  corresponding  values  of  p  ,  T0,  and  s0/R 
are  found  by  applying  Eqs.  (1-2),  (1-3),  and  (1-4)  in  Appendix  I.  These 
equations  are  shown  graphically  in  Figs.  1,  2,  and  3. 

In  the  tables  of  Ref.  5,  entries  are  clearly  indicated  by  dotted  lines 
for  which  dissociative  effects  exceed  1  percent.  On  this  basis,  the 
ranges  of  reservoir  conditions  are  intentionally  limited  to  exclude  such 
effects  almost  entirely. 


3.0  THROAT  CONDITIONS 


The  present  method  of  determining  the  throat  conditions  is  a 
generalization  of  the  one  developed  in  R,ef,  1.  In  writing  empirical 
equations  for  the  product  pu,  which  has  a  local  maximum  in  the  throat, 
it  is  obvious  to  start  with  the  conditions  of  constant  entropy  and  to  apply 
the  continuity  equation  for  the  conservation  of  energy,  which  is  essen¬ 
tially  a  relation  between  enthalpy  and  flow  velocity.  The  product  pu  can 
be  expressed  as  a  function  of  entropy  and  density  which,  at  constant 
entropy,  can  be  differentiated  with  respect  to  density. 

The  argument  begins  by  forming  the  function 

p2 u  =  2p2  (hD  --  h) 

Its  derivative  is  then 

=  4/j  ( h o  -  h)  -  2p2(dh/dp)s 
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This  derivative  vanishes  in  the  throat,  so  that 
2  (h0  -  h*)  -  P*  (dh/dp)s  =  0 

At  this  point  it  is  advantageous  to  introduce  the  identity 

WSp).  -  <W  (4S),  . 

On  substituting  Eq.  (5)  into  Eq.  (4),  there  follows 


/  d  log  h  \  * 

2  \d\ogp  j s 


(4) 

(5) 

(6) 


Equation  (6)  expresses  the  equilibrium  throat  condition  which  must  hold 
for  real  as  well  as  ideal  gases  and  its  evaluation  is  not  difficult. 

By  observing  the  identity 


and  that 


(  d  l°g  p  \ 

V  3  1osp/s 


=  yE 


where  k  is  the  dimensionless  enthalpy  function  hp/p  and  yg  is  the  isen- 
tropic  expansion  coefficient,  Eq.  (6)  becomes 

ho  =  h*  [(2m*  +  1  )/2m*  ]  (?) 

in  which  m  is  the  ratio  k/yp;  and  m*  is  its  value  in  the  throat.  Forget¬ 
ting  the  asterisks  for  a  moment,  the  function  comprising  the  entire 
right-hand  side  of  Eq.  (7)  can  be  tabulated  (or  plotted)  isentropically  as 
a.  function  of  h,  and  this  is  the  same  as  a  table  (or  plot)  of  h*  as  a  func¬ 
tion  of  h0  at  constant  entropy.  In  like  manner,  an  isentropic  table  (or 
plot)  of  the  right-hand  side  of  Eq.  (7)  as  a  function  of  density  yields  the 
relation  between  p*  and  h0.  These  relations  are  shown  graphically  in 
Figs,  4  and  5,  and  the  corresponding  empirical  equations  are  given  as 
Eqs.  (1-5)  and  (1-6)  in  Appendix  I.  The  numerical  value  of  the»product  , 
p*u*  is  then  given  by 

P*u*  =  P*  vm7-  h*)  (8) 
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Another  derivation  of  Eq.  (7)  is  given- in  section  4.  2.  .  It  should  also  be 
noted  that  the  present  parameter  m(  =  k /jj?)  is  akin  to  the  reciprocal  of  the 
parameter  m  used  in  Ref.  1. 


4.0  THE  SHOCK  CROSSING 


4.1  BASIC  PARAMETERS  AND  THEIR  INTERRELATION 

Although  the  procedure  being  set  forth  here  differs  somewhat  from 
the  one  developed  in  Ref.  1,  the  shock  crossing  again  involves  the  use 
of  four  dimensionless  parameters,  namely 

n  =  (h0  “  M/ho  =  1  h./ho  (9) 


r  =  pi/  Pi  =  Uj/ii! 

(10) 

hi  =  hiPi/pi 

(11) 

k2  =  h2p2/p2 

(12) 

In  the  present  instance,  the  fluid  in  front  of  the  shock  can  be  taken  as  a 
perfect  gas,  so  that  the  parameter  kj  is  substantially  constant,  its 
average  numerical  value  being  3,  4985  with  a  maximum  error  of  0.25  per 
cent.  Needless  to  say,  it  is  also  necessary  to  observe  the  three  con¬ 
tinuity  equations 

Pi  u,  =  p2U2  (13) 

Pi  + 

PiU*  -  P2  +  P2U2 

(14) 

2 

Ui  - 

t-  2H,  =  u/  +  2h2 

(15) 

These  four  definitions  and  the  continuity  equations  give  rise  to  the 
following  relations:  From  Eqs.  (9),  (10),  (1'3),  and  (15) 

UjJ  =  2nh0 

(16) 

u2  =  ru. 

(17) 

u/  =  2nr2h0 

(18) 

from  Eqs.  (9),  (15), 

and  (18) 

h,  =  (1  -  n)h0 

(19) 

h2  =  (1  -  nr2)h0 

(20) 
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from  Eqs.  (11)  and  (19) 

Pi  =  Pi/  k2  =  (1  -  n  )h0/5i  /kt  (21) 

and  from  Eqs.  (10),  (14),  and  (21) 

p2  =  Pi  +  piu/  (1  -  r)  (22) 

Pa/pi  =  1  +  2  nkj  ( 1  -  r )  /  ( 1  -  n)  (23a) 

From  the  definitions  of  kj  and  k.2,  the  pressure  ratio  across  the  shock 
can  also  be  written  in  the  form 

P2/P1  =  [kI  (1  -  nrs)j/[rk2(l  -  n  )J  (23b) 

The  following  four  relations  are  sometimes  useful  and  can  be  written  by 
inspection 

Pi/pi  =  (1  -  n)  (h0/ki) 

Pi/ Pa  =  r ( 1  -  n)  (ho/kt) 

Pa/Pi  =  [l  ~  n  +  2  nkx  ( 1  -  r )]  (ho/kj 
P2/Pa  =  r[l  ~  n  +  2 nkt  ( 1  -  r)](h0/k1) 

Finally,  by  substituting  Eqs,  (10),  (16),  (18),  (19),  and  (20)  into 
Eq.  (23),  there  follows 

(1  -  nr2)ki  —  2nr(l  -  r)kik2  —  r(l  -  n)k2  t=  0  (24) 

which  interrelates  the  four  dimensionless  parameters  across  the  shock 
and  can  easily  be  shown  to  be  an  alternate  form  of  the  Hugoniot  equation. 
This  form  of  the  Hugoniot  equation  is  especially  useful  in  the  present 
application  because  ki  is  substantially  constant  and  because  k2  can  be 
determined  from  observations  at  the  stagnation  point  downstream.  By 
rearranging  terms,  Eq.  (24)  is  a  quadratic  in  r,  namely: 

nkj  (2k2  -  l)r2  -  k2  (2nk,  +  1  —  n)r  +  ki  =  0  (25) 

The  Hugoniot  equation  is  often  written  in  the  more  familiar  form 

2  h2  -  2  h2  =  ( Pa/Pi )  +  (p  t/Pt)  ~  (pi/Pi )  “  (pi/p2) 

On  substituting  from  Eqs.  (11)  and  (12)  and  from  the  four  equations 
inserted  between  Eqs.  (23b)  and  (24),  this  form  of  Hugoniot's  equation 
becomes  a  simple,  but  rather  long,  expression  in  which  12  terms  annul 
one  another  in  pairs.  Equation  (25)  then  follows  by  assembling  the 
coefficients  of  r  and  r2. 
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For 'a  perfect  gas,  ki  =  k2  =  k,  and  with  this  restriction  Eq,  (24) 
reduces  itself  to  nr  (2k  -  1)  =  1.  Finally,  when  k  is  given  the  numerical 
value  3.5,  Hugoniot's  equation  for  a  perfect  gas  is  simply  6nr  =  1. 


4.2  THE  FREE-STREAM  MACH  NUMBER 

The  Mach  number  in  the  free -stream  is  a  function  of  the  parameters 
n  and  mi  (  =  ki/ygi).  The  free-stream  velocity  is  given  by 

u/  =  2nh0 

and  the  square  of  the  speed  of  sound  is 

ai  =  Y Pi/ Pi  ~  (1  ~  n ) h o Xjj i  / ki 

The  ratio  of  these  two  equations  is 

M,2  =  2  nm,/ (1  -  n ) 
or 

n  =  M1V(M12  +  2  m,) 

If  the  medium  can  be  construed  to  be  a  perfect  gas  for  which  k]_  =  3.  5 
and  YE1  =7=1.4,  so  that  mi  =  2.  5,  then 

M,2  =  5n/(l  -  n) 

Such  relations  are  also  valid  in  the  throat  where  the  Mach  number 
is  equal  to  unity,  so  that 

n*  =  1/  ( 2  m*  +  1) 


and  since 


it  follows  that 


1  -  n*  =  2m*/ (2m*  +  1) 


h*  =  h0(l  -  n*) 


h0  =  h*  (2m*  +  1 ) / 2 m* 


as  is  set  forth  in  Eq.  (7). 


4.3  IMPLICATIONS  OF  ISENTRQPIC  FLOW  BEHIND  NORMAL  SHOCK 

In  the  hypervelocity  regime,  the  thermodynamic  parameters  behind 
the  normal  shock  change  very  little  in  going  from  the  back  of  the  shock 
to  the  stagnation  point.  In  a  typical  run  at,  say,  4000°K  and  at  Mach  12, 
the  enthalpy  increases  by  about  1  percent,  the  density  by  about  5  percent, 
and  the  pressure  by  about  6  percent.  The  change  in  the  dimensionless 


8 


AEDC-TDR-64-50 


enthalpy  function  from  k2  to  k0'  is  hardly  perceptible.  Under  such  con¬ 
ditions  the  ratio  of  the  specific  heats  is  substantially  constant. 


This  renders  it  permissible  to  regard  the  fluid  behind  the  shock  as 
a  perfect  gas  of  some  kind,  especially  since  the  numerical  value  of  the 
specific  heat  ratio  is  immaterial.  Following  a  procedure  commonly 
adopted  in  the  literature  (for  example.  Ref.  10)  and  as  is  carried  through 
in  detail  in  Ref.  1,  the  argument  rests  on  the  isentropic  flow  equation 
for  a  compressible  gas 


Po  =  p2  1  + 


M, 


y/(y~  i ) 


(26) 


The  right-hand  side  of  this  equation  is  expansible  in  a  rapidly  con¬ 
vergent  series.  By  discarding  negligible  terms  and  applying  the  rela¬ 
tion  for  a  compressible  gas,  (y/2)M2^  =  ( 1 /2)p2U22/P2j  the  specific  heat 
ratio  drops  out  entirely,  and  Eq.  (26)  then  becomes 

Po  =  p2  +  (1/2)  p2u22  (1  +  w)  (27) 

p 

where  w  is  an  abbreviation  for  the  quantity  (1/4)M2  .  Then,  by  applying 
the  mass -flow  and  momentum  continuity  equations,  Eq.  (27)  becomes 

Po  =  Pl  +  ( 1/2 )  pi  u/  [2  -  r  ( 1  -  w)] 

The  density  ratio  is  always  less  than  0.  16  and  the  numerical  value 
of  (1-w)  is  slightly  less  than  unity.  It  follows  that  the  product  r(l-w) 
can  be  more  than  10  percent  in  error  before  its  effect  on  the  quantity 
in  brackets  exceeds  1  percent.  This  is  an  important  fact,  and  as  will 
presently  be  seen,  good  use  can  be  made  of  it. 

The  quantity  w  is  small  and  it  always  appears  in  the  expression 
[2  -  r(l-w)J  or  its  equivalent.  Since  the  density  ratio  is  also  small,  it 
follows  that  any  reasonably  correct  approximation  for  the  value  of  w  is 
adequate.  Through  the  range  of  conditions  for  which  this  program  is 
intended,  the  Mach  number  behind  the  shock  will  be  in  the  range  from 
0.  28  to  0.  38.  The  corresponding  values  of  w  lie  in  the  range  from  0.02 
to  0.  04,  so  that  0.  03  may  be  taken  as  an  average.  On  combining 
Eqs.  (16),  (21),  and  (27),  there  follows 

[  (  1  —  n  )  +  nkj  (2  —  0.97r)]ho  (28) 

in  which  the  enthalpy  is  expressed  in  atm/amagat  units.  The  parameter 
(1-n)  increases  from  about  0.  01  to  0.  05  as  the  free -stream  Mach  num¬ 
ber  decreased  from  20  to  10.  Inasmuch  as  this  program  is  intended  for 
use  at  Mach  numbers  from  10  upwards,  Eq.  (28)  can  be  viewed  from  the 
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point  of  view  that  the  numerical  value  of  the  term  (1-n)  in  the  denomi¬ 
nator  does  not  exceed  0.  05.  The  remainder  of  the  expression  inside 
the  square  brackets  keeps  a  numerical  value  about  6.5.  If,  in  the 
extreme  case  at  Mach  10,  the  term  (1-n)  is  neglected,  the  resulting 
error  is  about  0.  8  percent,  and  it  decreases  as  the  Mach  number  in¬ 
creases.  If  for  the  moment  the  term  (1-n)  is  neglected,  Eq.  (28) 
becomes 

Pi  =  p'J  [nh0  (2  -  0.97  r)]  (29) 

The  omission  of  the  term  (1-n)  in  this  equation  is  adequately  taken  into 
account  in  the  ultimate  determination  of  the  free -stream  density. 


4.4  DETERMINATION  OF  THE  PARAMETERS  n  AND  r 


It  is  well  known  that  shock  crossing  calculations  for  real  gases 
involve  an  iteration  of  one  kind  or  another.  The  usual  procedure  is  to 
assume  a  numerical  value  of  some  thermodynamic  or  flow  parameter 
and  with  recourse  to  the  continuity  equations,  to  test  the  implications 
of  the  assumption  with  reference  to  known  stagnation  conditions.  Various 
assumptions  are  tried  until  a  set  of  results  are  found  which  are  mutually 
consistent  within  prescribed  limits. 


In  principle,  the  present  method  also  invokes  iteration,  but  the 
iteration  is  hidden  by  the  device  of  numerical  estimates  in  situations 
which  are  insensitive  to  the  accuracy  of  the  estimates.  Besides,  in¬ 
stead  of  reserving  the  stagnation  conditions  for  purposes  of  comparison 
at  the  end  of  an  iterative  cycle,  they  are  put  into  the  numerical  works 
at  the  beginning.  The  method  is,  based  on  the  simultaneous  solution  of 
three  equations  in  logarithmic  form,  as  follows:  In  front  of  the  shock, 
the  perfect  gas  law  can  be  written  in  the  form 

RT,  =  (1  -  n)h0/k1 
or 

log  T,  =  log  (1  -  n)  +  log  li0  -  log  ki  -  log  R  (30) 


in  which  hQ  is  expressed  in  atm/amagat  units.  Likewise,  Eq.  (29)  can 
be  written  in  the  form 

log  Pi  =  Lg  Po  “  log  n  —  log  h0  -  log  (2  -  0.97  r)  (31) 

Next,  to  use  the  fact  that  the  upstream  entropy  is  known,  it  is  advan¬ 
tageous  to  resort  to  an  empirical  equation  in  the  form  of  the  theoretical 
expression  for  the  entropy  of  a  perfect  gas,  namely, 

log  T\  =  0.17364  (s0/R)  +  0.39971  log  Pl  -  1.50950  (  32) 


10 


AEDC-TDR-64-50 


Using  the  numerical  values  of  R  and  k]_  and  eliminating  Ti  and  from 
these  three  equations,  there  follows 

log(l  -  n)  =  0.17364  (s0/R)  -  1.39971  log  h0 
+  0.39971  log  Po  -  0.39971  log  n 

-  0.39971  log  (2  -  0.97  r )  -  3.40202  (33) 

The  purpose  of  this  equation  is  to  yield  the  value  of  the  parameter  n 
in  terms  of  a  somewhat  inaccurate  value  of  the  term  (1-n).  If,  in  the 
relation 

n  =  1  -  (1  -  n)  (34) 

the  numerical  value  of  (1-n)  is  approximately  0.  05,  then  a  20  percent 
error  in  (1-n)  causes  only  a  1  percent  error  in  n.  It  should  be  noted 
that  the  right-hand  side  of  Eq.  (33)  contains  a  term  log  n  and  another 
involving  the  density  ratio.  If  0.  97  is  taken  as  an  average  value  of  n, 
it  can  never  be  more  than  2  percent  in  error  in  the  Mach  range  from  10 
to  22,  However,  since  this  term  has  the  coefficient  0.  39971,  the  use  of 
the  average  value  of  n  at  this  point  contributes  less  than  a  1  percent 
error  in  (1-n). 

The  density  ratio  across  the  shock  plays  only  a  minor  part  in  deter¬ 
mining  the  thermodynamic  state  upstream  of  the  shock,  this  being  the 
small  fox~ward  effect  of  the  compressibility  of  the  gas.  In  the  quantity 
(2  -  0.  97  r),  the  constant  0.  97  is  valid  to  1  percent.  On  the  other  hand, 
the  density  ratio  is  relatively  small,  that  is,  0.  16  or  less.  Conse¬ 
quently,  the  numerical  value  of  the  density  ratio  used  at  this  point 
needs  only  to  be  approximately  correct,  and  such  a  value  is  easily 
found.  Equation  (24)  can  be  rearranged  into  the  form 

k2r  ^2ukj  (1  —  r)  +  (1  —  n)J  =  k2  d  —  nr) 

If,  on  the  left-hand  side  of  this  equation,  the  term  (1-n)  is  assumed  to 
be  negligible  with  respect  to  the  term  2nki(l-r),  the  error  introduced 
at  Mach  10  is  less  than  1  percent.  With  this  change,  the  above  equa¬ 
tion  becomes 

2k2nr  (1  -  r)  =  1  —  nr2 

By  long  division,  the  ratio  (l-nr2)/(l-r)  generates  the  expression 

1  +r  +r2(l  -  n )  /  ( 1  —  r) 

in  which  the  third  term  is  clearly  negligible.  The  above  equation  is 
thus  further  reduced  to 

2k2nr  =  1  +  r  (35) 
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By  giving  the  parameter  n  its  average  value  at  this  point  and  using  the 
condition  that  lt2  is  very  nearly  equal  to  ko',  there  follows 

7  =  l/(  1.94  k0'  -  1)  (36) 

where  the  tilde  in  r  denotes  an  approximate  value  of  r.  The  dimension¬ 
less  quantity  k0'  is  determined  by  an  empirical  equation  in  terms  of 
enthalpy  and  pressure  as  independent  variables.  See  Eq.  (1-8)  in 
Appendix  I  and  Fig.  6.  Equation  (33)  now  attains  its  working  form 

log  N  =  0.17354  (s0/R)  -  1.39971  log  h0  +  0.39971  log  Po 

-  0.39971  log  (2  -  0.97  7)  -  3.39673  (37) 

where  the  symbol  N  is  used  briefly  for  the  slightly  inaccurate  value  of 
(1-n).  The  final  value  of  n  is  given  by 

n  =  1  -  N  (38) 

and  the  final  value  of  the  density  ratio  is  found  by  solving  the  quadratic 
Eq.  (25)  which  is  reproduced  as  Eq.  (1-12)  in  Appendix  I. 


5.0  BEFORE-SHOCK,  AFTER-SHOCK,  AND  STAGNATION  CONDITIONS 


The  thermodynamic  state  in  front  of  the  shock  depends  very  much 
on  the  parameter  (1-n)  which  is  easily  found  after  the  density  and  tem¬ 
perature  are  determined.  The  density  is  given  by  a  slight  modification 
of  Eq,  (28),  namely 


Pi 


_ kl  Po _ 

N  +  nkt  (2  -  0.97  r)  h0 


(39) 


in  which  the  substitution  of  N  for  (1-n)  eliminates  practically  all  the 
error  introduced  in  going  from  Eq.  (28)  to  Eq.  (29).  The  temperature 
may  now  be  found  from  Eq.  (32)  (see  Fig.  7),  and  the  final  value  of 
(1-n)  is  given  by 

(1  -  n)  =  RT.k./ho  <40 


The  combination  of  terms  in  Eq,  (33)  clearly  shows  what  is  involved 
in  attaining  a  specified  degree  of  accuracy  in  determining  the  thermo¬ 
dynamic  state  in  front  of  the  shock.  The  heart  of  the  problem  is  seen 
in  Eq.  (32)  in  which  the  temperature  is  expressed  in  terms  of  the  up¬ 
stream  entropy.  Let  it  be  supposed  that  the  density  is  known  without 
error.  The  quantity  log  T  p  then  depends  directly  on  the  product 
(0.  17364)(s0/R).  If  the  temperature  is  to  be  determined  to  1  percent, 
log  Tp  must  be  known  to  within  ±0.  004.  In  the  hotshot  tunnels  for  which 
this  program  is  intended,  a  value  of  s0/R  equal  to  30.  0  is  fairly  typical. 
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and  the  product  (0.  17364)(30,  0)  is  equal  to  5.  209.  It  follows  that  the 
quantity  5.  209  cannot  be  in  error  by  more  than  0.  004  and  that  the 
entropy  must  be  known  to  0.  08  percent.  For  this  reason,  special  care 
was  used  in  constructing  Eq.  (1-4)  in  Appendix  I,  in  which  entropy  is 
expressed  as  a  function  of  pressure  and  enthalpy.  The  average  devia¬ 
tion  of  Eq.  (1-4)  from  125  equitably  distributed  entries  in  the  tables  is 
0.  06  percent.  Only  20  of  these  entries  deviate  from  the  equation  by 
more  than  0.  1  percent,  and  all  of  them  lie  in  the  lower-most  quartile 
of  the  enthalpy  range. 

The  remaining  quantities  of  interest  in  front  of  the  shock  now  follow 
directly.  The  pressure  p^  is  given  by  the  perfect  gas  law.  The 
enthalpy  Hi  and  the  flow  velocity  ui  follow  from  Eqs.  (1-17)  and  (1-18) 
and  the  speed  of  sound  from  Eq.  (1-19)  (see  Fig.  8)  in  Appendix  I.  The 
Mach  number  is  the  ratio  ui/ai,  and  the  effective  area  ratio  is  equal 
to  p*u*/p!Ui.  The  viscosity  is  given  by  Eq.  (1-22)  (see  Fig.  9)  and  the 
Reynolds  number  by  Eq.  (1-23)  in  Appendix  I. 

Behind  the  shock,  the  density,  pressure,  flow  velocity,  and  enthalpy 
follow  from  Eqs.  (10),  (23),  (17),  and  (20)  in  the  text  or  from  Eqs.  (1-24) 
through  (1-27)  in  Appendix  I.  The  temperature  and  speed  of  sound  are 
given  by  Eqs.  (1-28)  and  (1-29),  the  Mach  number  M2  being  thereby 
determined. 

At  the  stagnation  point,  the  density  p0'  is  determined  because 
k0',  hD,  and  pG'  are  known.  Temperature,  entropy,  and  viscosity  are 
computed  from  Eqs.  (1-32)  through  (1-34);  the  equation  for  the  entropy 
is  shown  graphically  in  Fig.  10.  The  Fay-Riddell  heat-transfer  rate 
follows  by  applying  Eq.  (1-35)  in  Appendix  I.  To  reconcile  the  meas¬ 
ured  and  computed  rates  of  heat -transfer,  the  computing  machine  is 
instructed  to  vary  the  overall  enthalpy  until  the  desired  degree  of  agree¬ 
ment  between  them  is  attained.  The  version  of  this  program  in  use  in 
the'  von  Karman  Gas  Dynamics  Facility  (VKF)  specifies  agreement 
within  1  percent,  and  this  result  is  usually  achieved  with  only  two  itera¬ 
tions  on  the  computing  machine. 


6.0  CONCLUDING  REMARKS 


The  pressures  in  the  reservoir  and  at  the  stagnation  point,  as  well 
as  the  heat-transfer  rate  at  the  stagnation  point  downstream,  are 
recorded  electronically  as  functions  of  time.  The  duration  of  a  run  is 
of  the  order  of  50  milliseconds.  In  selecting  corresponding  readings, 
it  is  necessary  to  take  into  account  the  time  required' for  a  pressure 
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signal  to  travel  from  the  throat  to  the  test  section.  This  interval  of 
time  depends  on  the  average  speed  v  with  which  the  signal  travels  and 
the  distance  from  the  throat  to  the  test  section.  This  average  speed 
can  be  written  in  the  form 

v  =  K  u* 

in  which  K  is  a  tunnel  constant  which  depends  on  the  expansion  angle  of 
the  nozzle,  the  throat  diameter,  and  the  thickness  of  the  boundary  layer. 
The  signal  is  assumed  to  be  progressing  through  successive  elements 
of  volume  in  the  tunnel  defined  by  parallel  planes  normal  to  the  longi¬ 
tudinal  axis.  In  any  such  element  of  volume,  the  velocity  of  the  signal 
is  given  by 

u,'  +  a;  =  uj  (1  +  1/Mi )  =  u*  (ui/u*)  (1  +  1/Mi ) 

where  the  subscript  identifies  the  successive  elements  of  volume 
serially.  The  constant  K  is  then  the  average  value  of  the  quantity 
(ui/u*)(l  +  1/Mi)  which,  in  fact,  is  very  nearly  constant  along  the  length 
of  the  tunnel  except  for  a  very  short  distance  at  the  throat.  Besides, 
the  numerical  value  of  K  is  surprisingly  insensitive  to  the  presumed 
thickness  of  the  boundary  layer.  In  the  case  of  Tunnel  H  (formerly 
known  as  Hotshot  2)  in  the  VKF  with  an  assumed  boundary-layer  thick¬ 
ness  of  half  the  radius  of  the  test  section,  the  empirical  expression 
for  K  is 

K  =  2.55  +  0.20  d* 

which,  in  practice,  needs  to  be  evaluated  only  to  two  significant  figures. 
The  quantity  d*  is  the  diameter  of  the  throat  in  inches. 


REFERENCES 


1.  Grabau,  M. ,  Humphrey,  R.  L.  and  Little,  W.  J.  "Determination 

of  Test-Section,  After-Shock,  and  Stagnation  Conditions  in  Hot- 
shot  Tunnels  Using  Real  Nitrogen  at  Temperatures  from  3000 
to  4000°K.  "  AEDC-TN-61-82,  July  1961. 

2.  Lukasiewicz,  J.  ,  Jackson,  R.  ,'  van  der  Bliek,  J.  A.  ,  Harris,  W.  G.  , 

and  Miller,  R,  M.  "Development  of  Capacitance  and  Inductance 
Driven  Hotshot  Tunnels.  "  AEDC-TN-60-222,  January  1961. 

3.  Lukasiewicz,  J. ,  Jackson,  R.  and  Whitfield,  J.  D.  "Status  of 

Development  of  Hotshot  Tunnels  at  the  AEDC.  "  Presented  at 
the  AGARD  Meeting  on  "High  Temperature  Aspects  of  Hyper¬ 
sonic  Flow,  "  Training  Center  for  Experimental  Aerodynamics, 
Rhode-Saint-Genese,  Belgium,  April  1962. 


14 


AEDC-TDR-64-50 


4.  Fay,  J.  A.  and  Riddell,  F.  R.  "Theory  of  Stagnation  Point  Heat 

Transfer  in  Dissociated  Air.  "  Journal  of  the  Aeronautical 
Sciences,  Vol.  25,  No.  2,  pp.  73-85,  1958. 

5.  Smith,  C,  E.  ,  Jr.  "Thermodynamic  Properties  of  Nitrogen.  " 

Lockheed  Technical  Report  6-90-62-111,  December  1962. 

(cf.  Foreword) 

6.  Hilsenrath,  J.  and  Klein,  M.  "Tables  of  Thermodynamic  Proper¬ 

ties  of  Nitrogen  in  Chemical  Equilibrium,  Including  Second 
Virial  Corrections  from  2000  to  15,  000°K.  "  AEDC-TDR- 
63-162.  March  1964. 

7.  Little,  W.  J.  and  Neel,  C.  A.  "Tables  of  the  Thermodynamic 

Properties  of  Nitrogen  from  100  to  1500°K.  "  AEDC-TDR- 
62-170,  September  1962. 

8.  Amdur,  I.  and  Mason,  E.  A.  "Properties  of  Gases  at  Very  High 

Temperatures."  The  Physics  of  Fluids,  Vol.  1,  No.  5, 
pp.  370-383,  1958. 

9.  Lewis,  C.  H.  and  Burgess,  E.  G.  ,  III.  "Empirical  Equations  for 

the  Thermodynamic  Properties  of  Air  and  Nitrogen  to 
15,  000°K."  AEDC-TDR-63-138,  July  1963. 

10.  Liepmann,  H.  W.  and  Roshko,  A.  Elements  of  Gasdynamics, 

John  Wiley  &Sons,  1957. 


15 


.4 


7.8 


8.0 


7.6 

log  'KQ 

i  Entropy  in  the  Reservoir  as  a  Function  of 
i  Reservoir  Pressure  and  Overall  Enthalpy 


AEDC-TDR-64-50 


APPENDIX  I 

DETAILS  OF  COMPUTATIONAL  PROCEDURE 

Note  the  relations 

h  =  8.722  (10)5  h 
u  =  9.339  (10)J  a 

Where  h  is  expressed  in  ft2 /sec2  and  u  in  ft/sec. 

Calculate  h0  from  the  measured  values  of  q,  Ra  and  Po  . 

h0  -  1.459  (10)5  (q  Vld/VTT)  +  0.7750  (10)7  (1-1  ) 

Calculate  p0  from  h0  and  p0 

logp0  =  2.2415672  (10)1  -  7.4788666  ( log  h0 ) 

-  1.5430265  do)"1  (log  p0)  +  8.5032627  (10)"1  (logh0)2 

+  4.7065564  (10)“'  (logh0)  (logp0)  ~  3.5687463  (10)"1  (logPo)2 

-  3.5600516  (10ra  (log  ho)*  -  4.5287166  (10)“ 2  (logh0)2  (log  Po) 

+  6.0866340  (10)"2  (log  h0)  (log  pof  -  1.9536922  (10)“ 2  (log  p0)3  (1-2) 

Calculate  T0  from  h0  and  p0 

To  =  (10)3  [r  1.6784483  (10 )3  +  7.0663450  (10 )2  (log|h0) 

-  1.2484939  (log  po)  -  9.9473139  ( 10 )‘  ( log  h0 )2 

+  2.3547617  (10)"’  (logh0)  (log  Po)  +  2.4996545  (10)"1  (log  Po)2 
+  4.6856478  (log  h0)3  -  1.6175003  (10)“ 2  (log  h0)2  (log  Po) 

-  2.9673076  (lOT3  (log  h0)  (log^po)2  -  4.3807464  (10)" 2  (log  Po)3  ]  .  ( 1-3  ) 

Calculate  s0/R  from  h0  and  Po 

s0/R  =  (10)1  [  1.2995109  -  3.8889476  ( 10)" 1  ( log  h0 ) 

-  3.9336637  (10)” 1  (log  Pb)  +  8.3586816  ( 1 0 )” 2  (logh0)2 

+  2.7878653  (10)"2  (log  h0)  (log  Po)  -  1.3109092  (10)"2  (logPo)2]  (1-4) 
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Calculate  h*  from  h0  and  s0/R 

log  h  *  =  ~  6.9866080  (10)1  +  1.5500832  (10 )‘  ( log  h  0 ) 

+  6.8621580  (10 )‘ (log  s0 /R)  -  2.6780798  ( log  h0 )2 
+  7.8302690  (log  h0)  (log  s0 /R)  -  6.7397050  (10)1  (log  s0/R)2 
„  +  1.6707170  (lOT1  (log  h0)3  -  8.5424000  do)-1  (log  h0 f  (log  s0 /R) 

+  2.2566900  (log  h0)  (log  s0 /R)s+  1.0606850  (10)2  (log  s0/R)3  (1-5) 

Calculate  p*  from  h0  and  s0/R 

logp*  =  9.0011196  (10)1  ~  2.1784020  (10)‘  (log  h0) 

-  8.9551180  (10)1  (log  s0 /R)  -  1.4084480  (10 )‘  ( log  h0 )2 

+  1.7752234  (10)2  (log  h0)  (log  s0/R)  -  3.8213446  (10)2  (log  s0/R)2 
4  6.6653368  (10)“'  (logh0)3  +  2.0314374.  (10)1  (logh0)2  (logs0/R) 

-  1.6678346  (10)2  (log  h0)  (log  s0 /R)2  +  3.6829143  (10 )2  (log  s0/R)3  (1-6) 


Calculate  p*  u*  from  h  *  and  p* 


(1-7) 


Calculate  k<(  from  h0  and  po 

ko  =  log  po  +  b2  (1-8) 

Where 

bt  =  -  0.0557  (log  h0/R)  +  0.2076 

1.1962  (log  h0  /R  -  4.1718) 

1  -  exp  [-  240  (log  K0  /R  -  4.1718)  ] 

b2  =  0.8405  (log  h0 /R)  +  0.6181 

3.2083  (  log  K0/R  -  4.2813  ) 

1  -  exp  [-  100  ( log  E0  /R  -  4.2813  )  ] 

2.4159  (log  K0  /R  -  4.4750) _ 

1  —  exp  [  —  80.2  ( log  Ii0  /R  —  4.4750  )  ] 
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Calculate  the  provisional  density  ratio  r  from  ko 
r  =  1  /  (1.94  k'  -  1) 

Calculate  the  quantity  N  from  s0  /R  ,  h0,  po  and  r 

log  N  =  0.17354  (s0/R)  -  0.23562  log  h0 

+  0.39971  log  p o  -  0.39971  log  (2  -  0.97  ¥)  -  3.39673 

Calculate  the  quantity  n  from  N 
n  =  1  -  N 


(1-9) 


(1-10) 


(1-11) 


Calculate  r  by  solving  the  quadratic 
k0'  ( 2  n  ki  +  1  -  n  ) 


2 

r  — 


n  k,  (  2  k0'  -  1  ) 


r  + 


1 


n  (  2  k0'  -  1  ) 


=  0 


(1-12) 


NOTE: 


The  desired  value  of  r  is  the  smaller  of  the  roots.  To  preserve  accuracy,  note  that  the 


smaller  of  the  roots  is  equal  to 


n  (  2  k0"  -  1 ) 


divided  by  the  larger  of  the  two  roots. 


Calculate  p1  from  po ,  N,  n,  r,  h0  and  kx 

p  =  - - 

N  +  n  k»  (2  -  0.97  r)  [E0  /8.722  (  1 0  )5  ] 

Calculate  Ti  from  s0/R  and  pt 

log  Tt  =  1.7364  (10)"1  (s0/R)  +  3.9971  (10)"1  (log  Pl)  -  1.5095 

Calculate  the  quantity  (1  -  n)  from  T,,  k,  and  h0 
1  -  n  =  8.722  (10)5  RT,  ks / ET0 

Calculate  pi  from  and  T! 

Pi  =  Pi  RT! 

Calculate  h,  from  (1  —  n)  and  h0 

h,  =  (1  -  n)h0 


(1-13) 


(1-14) 


(1-15) 


(1-16) 


(1-17) 
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Calculate  Ui  from  n  and  h0 

U;  =  V  2  n  li0  (1-18) 

Calculate  aj  from 

If  Tt  ?  400°K  * 

a\  =  6.6883  (10)1  V  T, 

If  Tt  ?  400X 

5,  =  1.1055  (10)3  [-  2.3537  (10)“*  +  6.4129(10)“*  JtT  -  1.2988  (10)“4  ij  (1-19) 


Calculate  Mt  from  at  and  Uj 

Mi  =  /  aj  (I-?0) 

Calculate  A/A*  from  p*  u*,  pt  and  it! 

A/A*  =  p*u*//>,u1  (1-21) 


Calculate  pt  from  Tj 


1(1.1172X10)-  _^~r(~!T)/,J 


(1-22) 


NOTE: 

When  T,  100X  ,h  =  4.62  x  10“”  rl\ 

Calculate  Re! /ft  from  pt  ,  u,  and  fa 

Ret  _  7.806  (  10)“*  p,  u, 

ft  “  Pi 

Calculate  p2  from  p2  and  r 

Pi  =  P:  /r 

Calculate  p2  from  pt  ,  n  ,  ki  and  r 

p2  —  Pi  [l  +  2  n  k]  (1  —  r )  /  ( 1  -  n)J 


(1-23) 


(1-24) 


(1-25) 
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Calculate  u2  from  r  and  u3 

u2  =  r  ?!  (1-26) 

Calculate  h2  from  n  ,  r  and  hQ 

h2  =  (1  -  n  r2)  h0  (1-27) 

Calculate  T2  from  p2  and  h2 

T2  =  (10)3  [  b3  log  p2  4  b4  ]  (1-28) 

Where 

b  _  1.4125  (log  h2/R  -  4.1587  ) 

,  1  —  exp  [—67  (log  h2/R  —  4.1587  )] 

-  1-0759  (log  h2/R  -  4.4062  ) 

1  ~j  exp  [-56  (log  h2/R  -  4.4062  )] 

b4  =|3.72456  (log  h2/R)  -  12.46890 

4-26469  (log  h2/R  -  4-04565) 

1  —  exp  [  —  28.1  ( log  h2  /R  —  4-04565  )  ] 

-5.41781  (log  h2/ft  -  4.37536  ) 

1  -  exp  [-33  (log  hj/R  -  4.37536  )] 

Calculate  a2  from  p2  and  p2 

a2  =  1.1055  (10)3  [3.1491  log  p2  4-  b,, 

0.4808  ( log  p2  -  c,  ) 

1  -  exp  [—100  ( log  p2  —  Cj )  ] 

—  1.2419  log  p2  —  c2 

1  4-  exp  [-  27.5  ( log  p2  -  cs  )  ] 

4-  0.0553  exp  [~  100  (log  p2  -  c4)  ]  ]  (1-29) 

Where 

b5  =  -3.1491  log  p2  -  0.1167  c3  =  1.02  log  p2  +  1.21 

cs  =  0.9917  log  p2  +  1.0003  c4  =  log  p2  +  0.7397 

c2  =  1.3697  log  p2  +  1.5383 
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Calculate  M2  from  u2  and  a2 
M2  =  u2/a2 

Calculate  p0'  from  k0',  Po  and  h0 

Po  =  k0'  [  8.722  (10)5  po  /  h0  ] 

Calculate  T0'  from  po  and  h0 
Use  Eq.  (1-28) 

Calculate  s0'/R  from  p0'  and  h0 
s07R  =  b6  log  p0'  +  b7 

Where 

b6  =  -  0.0231  log  h0/R  -  2.2089 

-  1.2157  (log  h0  /R  -  4.2818) _ 

1  -  exp  [  -40.5  (log  h0  /R  -  4.2818)] 

-  4.6987  (log  h0 /R  -  4.5818) _ 

1  -  exp  [  -  26.4  (log  h0  /R  -  4.5818)  ] 

b7  =  9.0085  log  h0/R  -  4.7282 

1.7300  (log  h0  /R  -  4.2056  ) _ 

+  1  -  exp  [  -173  (log  h0/R  -  4.2056  )} 

Calculate  p0'  from  T0' 

Po  =  1.1172  (10)-5  [  1.0256  +  1.4223  (10)~3  T0'  -  1.8136  (10)“'  (TV)* 
Calculate  the  quantity  q  VRa  from  p0' ,  p0',  h0  ,  p0'  and  p2 
qVlU  =  4.2519  (lO)"4  (p0')OU  ( po  )°'15  [h0  -  3.3469  (10)'] 

,  ,  \  0 .25  /  /.0.10 

(po  “  Pi)  (Po) 


(1-30) 


(1-31) 


(1-32) 


]  (1-33) 


(1-34) 
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APPENDIX  II 

COMPUTER  PRINTOUTS  FOR  TYPICAL  RUNS 
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NITROGEN  Q  INPUT 


RUN  NUMBER 


1 


1964  NITROGEN  0  OOT  PROGRAM 

REFERENCE  DENSITY  0.07606410  LBS/CUBIC  FOOT 
RESERVOIR  CONDITIONS 


T  SEC 

PO  PS  I 

RHO-O  ATM 

TO  OEG  K 

HO  SOFT/SGSEC 

SOR 

.0000 

25000 . 

141.85 

2628. 

.357164-06 

24.031 

.0000 

25000  . 

129.65 

2935. 

.40042+08 

24.512 

.0000 

25000. 

88.85 

4590. 

.63681+08 

26.504 

co 


RADIUS  .50  INCHES  DATA  PROGRAM  006A 


STAGNATION  CONDITIONS  BEHIND  NORMAL  SHOCK 


O-DOT 

PO+PSI 

RHO-O+ATM 

TO  +  DEG  K 

SO  +  /R 

200.00 

8.000 

.0535110 

2766. 

32.35 

200.00 

6.000 

.0361977 

3082. 

33.03 

200. CO 

2  .  Cv/Q 

.0083140 

4418. 

36.20 
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NITROGEN  0  INPUT  RUN  NUMBER  .  1  RADIUS  .50  INCHES  DATA  PROGRAM 

1964  NITROGEN  0  DOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LBS/CUBIC  FOOT 
FREE  STREAM  CONDITIONS  MODEL  LENGTH  12.000  INCHES 


T  SEC 

P-1' PSI 

RHO-1  ATM 

T— 1  DEG  K 

U-l  FT/SEC 

M—  1 

0-1  PSI 

RE/INCH 

RE/MODEL 

A/ A* 

PO+/PO 

.0000 

.0252838 

.0072945 

64.42 

6366. 

15.58 

4.298 

133377. 

1600522. 

5318. 

.3200-03 

.0000 

.0173892 

.0043666 

66.41 

86  66. 

16.27 

3.220 

91473. 

1097675. 

7177. 

.2400-03 

.0000 

.0042549 

.0010075 

78.50 

11207. 

18.91 

1.065 

20252. 

243026. 

22386. 

.8000-04 

00 

on 


006  A 


V-l 

1148-01 

1447-01 

3576-01 
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RUN  NUMBER 


NITROGEN  Q  INPUT 
1964  NITROGEN  Q  DOT  PROGRAM 

REFERENCE  DENSITY  G .  0*78064!  0  LBS/CUBIC  FOOT  . 


CONDITIONS  BEHIND 

NORMAL 

SHOCK 

T  SEC 

P-2  PSI 

RHO-2  ATM 

7-2  DEG 

-0000 

7.38404 

.050446 

2716. 

.0000 

5.54387 

•034141 

302  7. 

.0000 

1.86378 

-007879 

4373. 

H-2  SQFT/SQSEC 
-401100+02 
-449936+02 
•  718349+02 


RADIUS 


U-2  FT/SEC 

1209. 74 

1263.74 
1433.11 


.50 


M-2 
.3608 
.  3565 
.2314 


INCHES 


DATA  PROGRAM 


CO 

05 


006A 
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NITROGEN  Q  INRUT  RUN  NUMBER  2 

1964  NITROGEN  0  DOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LBS/CUBiC  FOOT 
RESERVOIR  CONDITIONS 


RADIUS 


.  50 


INCHES 


DATA  PROGRAM 


STAGNATION  CONDITIONS  BEHIND  NORMAL 


T  SEC 

PC  PS  I 

RHO-O  ATM 

TO  DEG  K 

HO  SOFT/ SL SEC 

SOR 

C-DOT 

PO+PSI 

RHU-Q+ATM 

TO+DEG  K 

SO  +  /R 

.0000 

3C00C  . 

182.56 

2305. 

.31521+Ct 

23.297 

170.00 

8.000 

.0599333 

2474. 

31.86 

.0000 

30000. 

167.86 

2566. 

.35198+08 

23-757 

170-00 

6.000 

.0406906 

2729. 

32.58 

.coco 

30C0C . 

117.66 

3984. 

. 55292+G« 

25.682 

170.00 

2.  COO 

.0091776 

4034. 

35.51 

.0000 

3DGC0. 

91.43 

5332. 

.74984+08 

27.017 

170.0C 

■  1.000 

.0038748 

4700. 

37.76 

CO 

-J 


006A 


SHOCK 
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NITROGEN  Q 

INPUT 

RUN 

NUMBER  2 

RADIUS 

.50 

INCHES 

DATA  PROGRAM  006A 

1964  NITROGEN  Q  DOT 

PROGRAM 

REFERENCE 

DENSITY  0. 

0730641 C 

LBS/CUbIC  FOOT 

FREE 

STREAM  CONDITIONS 

MODEL 

LENGTH  12.000  INCHES 

T  SEC 

P-1  PSI 

RHO-1  ATM 

T-l  DEG 

K  U-l  FT/SEC  M-l 

Q-l  PSI 

RE/INCH 

RE/MODEL 

A/A* 

PO+/PO 

V-l 

.0000 

.0224080 

.0082533 

5C  .46 

7e69.  16.56 

4.302  • 

181199. 

2174390. 

6173. 

.2667-03 

.1047-01 

.0000 

.0153770 

.0055293 

51.69 

8321.  17.30 

3.223 

125329. 

1503951  . 

8358. 

.2000-03 

.1315-01 

.0000 

.0037398 

.0011620 

59.82 

10452.  20.21 

1.069 

28587. 

343047. 

26421. 

.6667-04 

.3215-01 

.0000 

.0015519 

.0004234 

68.13 

12184.  22.07 

.529 

10660. 

127924. 

53979. 

.3333-04 

.5751-01 

OJ 

03 
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NITROGEN  Q 

INPUT 

RUN 

NUMBER  2 

1964  NITROGEN  0  DOT 

PROGRAM 

REFERENCE 

DENSITY  0 

.07806410 

LBS/CUBIC  FOOT 

CONDITIONS  BEHIND 

NORMAL 

SHOCK 

' 

T  SEC 

P-2  PSI 

RH0-2  AIM 

T-2  DEG 

K 

H-2  SOFT /SQSEC 

.0000 

7.37343 

.056445 

2432. 

.353806+02 

.0000 

5.53589 

. 033341 

2680. 

.395301+02 

.0000 

1.85612 

.008674 

3983. 

.622694+02 

.0000 

.93893 

.003688 

4671. 

.848496+02 

RADIUS  ,50  INCHES  DATA  PROGRAM  006A 


U-2  FT/SEC 

H-2 

115C.54 

.3623 

1200.01 

.3601 

1400.13 

.2680 

1398.69 

.1856 

Co 

CO 


AEDC-TDR-64-50 


NITROGEN 

0  INPUT 

RUN 

NUMBER 

3 

RADIUS 

.50 

INCHES 

DATA 

PROGRAM 

1964  NITROGEN  0 

DOT  PROGRAM 

' 

REFERENCE  DENSITY 

0.07306410 

LBS/CUB IC  FOOT 

RESERVOIR  CONDITIONS 

STAGNATION  CONDITIONS  BEHIND  NORMAL 

T  SEC 

PO  PS  I 

RHO-O  ATM 

TO  DEG  K 

HO  SQF T /SijSEC 

SOR 

O-DOT 

PO+PSI 

RHO-O+ATM 

TO+DEG  K 

S0  +  /R 

.0000 

20000  . 

139.60 

2135. 

.28427+08 

23.344 

130.00 

6.000 

.0493808 

2269. 

31.74 

.0000 

20000. 

93.35 

3250. 

.44105+08 

25.177 

130.00 

2 . 000 

.0110312 

3372. 

34.56 

.0000 

20000. 

77.05 

4306. 

. 59164+Ot, 

26.436 

130.00 

1.000 

.0044017 

4176. 

36.54 

.0000 

20000. 

70. 9" 

4724.' 

.65233+08 

26.860 

130.00 

.600 

. 0033385 

4395. 

37.26 

C06  A 


SHOCK 


rfs. 

O 


A  EDC-TD  R-64-50 


NITROGEN  0  INPUT  RUN  NUMBER  3  RADIUS  .50  INCHES  DATA  PROGRAM  006A 

1964  NITROGEN  0  DOT  PROGRAM 

REFERENCE  DENSITY  0.078064X0  LBS/CUB IC  FOOT  , 

FREE  STREAM  CONDITIONS  MODEL  LENGTH  12. COO  INCHES 


T  SEC 

P-1  PS  I 

RHO-1  ATM 

T-l  DEG  K 

U-l  FT/SEC 

M-l 

0-1  PSI 

RE/INCH 

RE/MODEL 

A/  A» 

PC+/PO 

V-  1 

.0000 

.0176898 

.0068760 

47.82 

7469. 

16.15 

3.229 

151223. 

1814672. 

5664  • 

.3000-03 

.1117-01 

.0000 

.0042203 

.0014635 

53.60 

9328. 

19.05 

1.C72 

35862. 

430344. 

17621  . 

.1000-03 

.2706-01 

.0000 

.0017358 

.0005414 

59.59 

10816. 

20.95 

.533 

13837. 

166033 . 

35942. 

.5000-04 

.4792-01 

.0000 

.0013056 

.0003913 

62.02 

11361. 

21.57 

.425 

10092. 

121104. 

45116. 

.4000-04 

.5777-01 

AE  DC-TDR-64-50 


NITROGEN  0  INPUT 


RUN  NUMBER 


3 


1964  NITROGEN  Q  DOT  PROGRAM 


REFERENCE  DENSITY  0.0-7836410  LBS/CUBIC  FOOT 
CONDITIONS  BEHIND  NORMAL  SHOCK 


T  SEC 

P-2  PS1 

RHO-2  ATM 

T-2  OEG  K 

H-2  SOFT/SOSEC 

.0000 

5.52479 

.046484 

2231. 

.318925+02 

.0000 

1.84942 

.010453 

3316. 

.495901+02 

.0000 

.93028 

.004166 

4131. 

-  6  67006  +02 

.0000 

.74714 

.003167 

4357. 

.736617+02 

CO 


■  RADIUS 


50 


INCHES 


DATA  PROGRAM  006A' 


U-2  FT/SEC  M-2 

1104.83  .3644 

1306.04  .3379 

1405.66  .2229 

1403.62  .1958 


AEDC-TDR-64-50 


NITROGEN  C  INPUT  RUN  NUMBER  4 

1964  NITROGEN  0  DOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LBS/CUBIC  FOOT 
RESERVOIR  CONDITIONS 


RADIUS 


.50 


INCHES 


DATA  PROGRAM 


STAGNATION  CONDITIONS  BEHIND  NORMAL 


T  SEC 

PO  PSI 

RHO-C  ATM 

TO  DEG  K. 

HO  SQFT/SCSEC 

SOR 

U-DOT 

PO+PSI 

RHQ-G+ATM 

TO+OEG  K 

S0  +  /R 

.0000 

4000  . 

33.16 

2139. 

.27050+02 

24.953 

50.00 

1 .000 

.0086272 

2178. 

33.34 

.0000 

4000. 

30.42 

2340. 

. 29859+08 

25.350 

50.00 

.800 

.CC 63165 

2363. 

33.95 

.0000 

4000. 

27.65 

2585. 

.33279+08 

25.709 

50.00 

.600 

. 0042936 

2594. 

34.66 

.0000 

4000. 

24.02 

2994.  . 

.39017+08 

26.440 

50.00 

.400  • 

.0024850 

3005. 

35.69 

.0000 

4000. 

18.56 

3912. 

. 5 1968+06 

27.633 

50.00 

.200 

.0009811 

3781. 

37.54 

CO 


006A 


SHOCK 


AEDC-TDR-64-50 


NITROGEN  Q  INPUT  RUN  NUMBER  4  RADIUS  .50  INCHES  DATA  PROGRAM  006A 

1964  NITROGEN  Q  DOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LSS/CU8 I C  FOOT 
FREE  STREAM  CONDITIONS  MODEL  LENGTH  12.000  INCHES 


T  SEC 

P-1  PSI 

RHO-1  ATM 

T-l  DEG  K 

U— 1  FT/SEC 

K— 1 

Q-l  PSI 

'RE/INCH 

RE/MODEL 

A/ A» 

PD+/P0 

V—  1 

.0000 

.0029380 

.0012045 

45.34 

7286. 

16.18 

.538 

27256. 

327067. 

6890. 

.2500-03 

.2637-01 

.0000- 

.0021873 

.000.8710 

46.68 

7660. 

16.76 

.430 

20126. 

241514. 

8631. 

.2000-03 

.3179-01 

.0000 

.0014920 

.0005846 

47.44 

8C93. 

17.5  7 

.322 

14042. 

168502. 

11556. 

.1500-03 

.3989-01 

.0000 

.0008740 

.0003313 

49.04 

8/71. 

18.73 

.215 

8342. 

100108. 

17440. 

.1000-03 

.5516-01 

.0000 

.0003539 

.0001235 

53.27 

10136. 

20.76 

.107 

3309. 

39708. 

35148. 

.5000-04 

.9712-01 

AEDC-TDR-64-50 


NITROGEN  0 

INPUT 

RUN 

NUMBER  4 

1964  NITROGEN  Q  DOT 

PROGRAM 

REFERENCE 

DENSITY  0 

.07806410 

IBS/CUBIC  FOOT 

CONDITIONS  BEHIND 

NORMAL 

SHOCK 

T  SEC 

P-2  PSI 

RHO-2  ATM 

T-2  DEG 

K 

H-2  SQFT/SQSEC 

.0000 

.92053 

.008120 

2141. 

.303442+02 

.0000 

.73718 

.005948 

2324. 

.335126+02 

.0000 

.55350 

.004050 

2549. 

.373730+02 

.0000 

.36959 

.002343 

2953. 

.438519+02 

.0000 

.  18569 

.000927 

3741. 

.585378+02 

RADIUS 


50 


INCHES 


OATA  PROGRAM  006A 


U-2  FT/SEC 
1080.83 
1121.70 
1168.19 
1240.11 
1349.73 


M-2 

.3650 

.3619 

.3582 

.3444 

.2225 


Cn 


AEDC-TDR-64-, 


NITROGEN  Q  INPUT 


RUN  NUMBER 


5 


1964  NITROGEN  Q  DOT  PROGRAM 


REFERENCE  OENSITY 
RESERVOIR  CONDITIONS 


T  SEC 

PO  PS  I 

RHO-O  ATM 

.0000 

6000. 

76.16 

.0000 

6000  . 

54.68 

.0000 

6000  . 

42.81 

.0000 

6000  . 

39.77 

.0000 

60CC  . 

36.05 

.0000 

6000. 

31.18 

05 


0.07806410  LBS/CUBIC  FOOT 


TO  DEG  K 

HO  SQFT/SCSEC 

SOR 

1302. 

. 15907+0F 

22.417 

1880. 

.23630+08 

23.967 

2442. 

.31480+08 

25.122 

2642. 

.34281+08 

25.470 

2935. 

.38385+08 

25.936 

3424. 

.45270+08 

26.622 

RADIUS  .50  INCHES  DATA  PROGRAM  006A 


STAGNATION  CONDITIONS  BEHIND  NORMAL  SHOCK 


Q-DOT 

PO+PSI 

RHO-O+ATM 

TO+DEG  K 

SO+/R 

60 .00 

6.000 

.0833787 

1300. 

29.47 

60.00 

2.000 

.0194727 

1946. 

32.11 

60.00 

1.000 

.0075258 

2471. 

33.93 

60.00 

.800 

.0055774 

2664. 

34.49 

60.00 

.600 

.0037796 

2959. 

35.22 

60.00 

.400 

.0021733 

3440. 

36.23 

AEDC-TDR-64-50 


4^ 

-a 


NITROGEN  0 

INPUT 

RUN 

NUMBER  5 

RADIUS 

.50 

INCHES 

DATA  PROGRAM 

1964  NITROGEN  C  OCT 

PROGRAM 

REFERENCE 

DENSITY  0 

.07306410 

LBS/CUBIC  FOOT 

FREE 

STREAM  CCiMDITICNS 

MODEL 

LENGTH  12.000  INCHES 

T  SEC 

P-1  PSl 

RHO-1  ATM 

T-l  OEG 

K  U-l  F  T /SEC  M— 1 

0-1  PSI 

RE/INCH 

RE/MODEL 

A/ A» 

PO+/P j 

.0000 

.0261673 

.0124949 

41.90 

5557.  12.84 

3.243 

233337. 

28 C 0045. 

1760. 

.1000-02  . 

.0000 

.0063419 

.0027t>54 

42.6? 

6305.  15.58 

1.078 

62164. 

745970. 

5211. 

.3333-03  . 

.  .0000 

.0025262 

.0010308 

45.59 

7870.  17.43 

.537 

25157. 

300688. 

1 04 P  5  . 

.1667-03  . 

.0000 

.0013821 

.0007558 

46.28 

rt217.  18.06 

.43  0 

16894. 

226733. 

13143. 

.1333-03  . 

.0000 

.0012894 

.0005051 

47.45 

6701.  18.39 

.  322 

1,3039. 

156473. 

17596. 

.1000-03  . 

.0000 

.0007602 

.0002846 

49.64 

9457.  20.07 

.214 

7634. 

91607. 

26530. 

.6667-04  . 

006A 


V—  1  - 

7149-02 
1662-01 
2962-01 
3535-01 
4450-01 
6179-01 


A  EDC-TDR-64-50 


co 


NITROGEN  0 

INPUT 

RUN 

NUMBER  5 

RADIUS 

.50 

1964  NITROGEN  0  DOT 

PROGRAM 

REFERENCE 

DENSITY  0 

.07806410 

LBS/CUBIC  FOOT 

CONDITIONS  BEHIND 

NORMAL 

SHOCK 

T  SEC 

P-2  PSI 

RHO-2  ATM 

T-2  DEG 

K 

H-2  SQFT/SQSEC 

U-2  FT.'CEC 

M-2 

.0000 

5.492‘16 

.078279 

1262. 

.177872+02 

887.04 

.0306 

.0000 

1.83830 

.018314 

1910. 

.264875+02 

1027.59 

.3701 

.0000 

.92186 

.00  7088 

2429. 

.353412+02 

1144.62 

.3605 

.0000 

.73814 

.005255 

‘  2618. 

.385027+02 

1181.91 

.3574 

.0000 

.55423 

.003563 

2907. 

.431372+02 

1233.37 

.3488 

.0000 

.37009 

.002050 

3388. 

.509153+02 

1312.62 

.30  36 

INCHES 


DATA  PROGRAM 


006A 


A  E  DC-T D  R-64-50 


NITROGEN  0  INPUT 


RUN  NUMBER 


6 


1964  NITROGEN  Q  DOT  PROGRAM 

REFERENCE  DENSITY  5. 07806410  LBS/CUBIC  FOOT 
RESERVOIR  CONDITIONS 


T  SEC 

PC  PS  I 

RHO-O  ATM 

TO  DEG  K 

HO  SQFT/SCSEC 

SOP. 

.C000 

7000. 

33.93 

1363. 

. 1 6758+CE 

22.443 

.0000 

7000 . 

75.73 

1532. 

.18983+03 

22.929 

.0000 

7006. 

64.15 

1840. 

.23176+0e 

23.718 

.0000 

7000. 

44.44 

2737. 

.35716+08 

25.470 

.0000 

7000. 

34.79 

3560. 

.47299+03 

26.639 

co 


RADIUS 


50 


INCHES 


DATA  PROGRAM  006A 


STAGNATION  CONDITIONS  BEHIND  NORMAL  SHOCK 


Q-UOT 

PO+PSI 

RHO-O+ATM 

TO+DEG  K 

SO  +  /R 

100.00 

14.000 

.1856524 

1335. 

23.83 

100. CO 

10.000 

.1184994 

1587. 

29.65 

ICO. 00 

6.000  • 

.0593945 

1913. 

30.94 

100.00 

2.000 

.0134135 

2766. 

33.73 

100.00 

1.000 

.0052155 

3574. 

35.53 

AEDC-TDR- 64-50 


NITROGEN  Q 

INPUT 

RUN 

NUMBER  6 

RADIUS 

-50 

INCHES 

DATA  PROGRAM  006A 

1964  NITROGEN  0  DOT 

PROGRAM 

* 

REFERENCE 

OENSITY  0 

.07806410 

L3S/CUBIC  FOOT 

FREE  STREAM  CONDITIONS 

MODEL 

LENGTH  12.000  INCHES 

T  SEC 

P-1  PSI 

RHO-l  ATM 

T— 1  DEG 

K  U— 1  FT/SEC 

M-l 

0-1  PSI 

RE/INCH 

RE/KODEL 

A/A» 

PO+/PO 

V-l 

.0000 

.0876415 

.0279491 

58.42 

5676. 

11.10 

7.581 

382398. 

4588777. 

873. 

.2000-02 

.4831-02 

.0000 

.0554959 

.0175208 

58.87 

6C5  5. 

11.80 

5.407  . 

253718. 

3044612. 

1216. 

.1429-02 

.6302-02 

.0000 

.0273375 

.0085444 

60.56 

6709. 

12.89 

3.237 

133280. 

1599363. 

2021. 

.8571-03 

.9499-02 

.0000 

.0064582 

.0018240 

65.81 

8364. 

15.42 

1.074 

32644. 

391724. 

6135. 

.2857-03 

.2296-01 

.0000 

.0026125 

.0006841 

70.98 

9644. 

17.12 

.536 

13089. 

157071. 

12397. 

.1429-03 

.4025-01 

AEDC-TDR-64-50 


NITROGEN  Q  INPUT  RUN  NUMBER  6  '  RADIUS  .50  INCHES  DATA  PROGRAM 

1964  NITROGEN  0  DOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LBS/CU8 IC  FOOT 
CONDITIONS  BEHIND  NORMAL  SHOCK 


T  SEC 

P-2  PSI 

RHO-2  ATM 

T— 2  DEG  K 

H-2  SQFT/SQSEC 

U-2  FT/SEC 

M-2 

.0000 

12.83189 

.174531 

1346. 

.187395+02 

909.01 

.1565 

.0000 

9.17673 

.111449 

1550. 

.212454+02 

951.87 

.3804 

.0000 

5.51680 

.055904 

1877. 

.259689+02 

1025.34 

.3728 

.0000 

1.84660 

.012648 

2717. 

.401148+02 

1206.27 

.3586 

.0000 

.92585 

.004924 

3521. 

.532003+02 

1339.87 

.3068 

006A 


AEDC-TDR-64-50 


NITROGEN  0  INPUT 


RUN  NUMBER 


7 


1964  NITROGEN  0  DOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LP.S/CUi3lC  FOOT 
RESERVOIR  CONDITIONS 


T  SEC 

PO  PS  I 

RHO-1'.  ATM 

TO  DEG  K- 

HO  SOFT / SC S EC 

SOR 

.0000 

12500. 

114.43 

1695. 

.21690+03 

22. '794 

.0000 

12500. 

103.86 

1897. 

.74486+08 

23.280 

.0000 

125GC  . 

87.27 

2317, 

. 30354+C8 

24.154 

.ooco 

12500. 

60.70 

3493. 

.46902+03 

25.967 

CO 


RADIUS  .  .50  INCHES  DATA  PROGRAM  006A 


STAGNATION  CONDITIONS  8EHIND  NORMAL  SHOCK 


Q-DOT 

PO+PSI 

RHO-O+ATM 

TO+DEG  K 

SO  +  /R 

140.09 

14.000 

. 1470502 

1804. 

29.84 

140.00 

.  1C. 000 

.0941474 

2005. 

30.65 

140.00 

6  .  OGO 

.0465334 

2396. 

32.00 

140.00 

2.000 

.0104859 

355S. 

34.80 

A  EDC-T  DR-64-50 


RUN  NUMBER 


7 


RADIUS 


50 


INCHES 


NITROGEN  Q  INPUT 
1964  NITROGEN  0  DOT  PROGRAM 


REFERENCE 

DENSITY  0 

.07806410 

REE 

STREAM  CONDITIONS 

MODEL 

T  SEC 

P-1  PSI 

RHO-1  ATM 

T— 1  DEG  K 

.0000 

.0693751 

.0213607 

6C  .37 

.0000 

.0441154 

.0134529 

60.9  5 

.0000 

.0224525 

.0064704 

64.50 

.0000 

.0053350 

.0013808 

71.82 

LBS/CUBIC  FOOT 
LENGTH  12.000  INCHES 


U— 1  - FT/SEC 

M—  1 

Q-l  PSI 

RE/INCH 

6484. 

12.48 

7.559 

323037. 

6900. 

13.22 

5.392 

214444. 

7699. 

14.33 

3.228 

108747. 

9602. 

16.94 

1.072 

25994. 

DATA  PROGRAM  006A 


RE/MODEL 

A/ A» 

PO+/PO 

r— 4 

1 

> 

3876443. 

1534. 

.1120-02 

.5906-02 

2573331 . 

2155. 

.8000-03 

.7678-02 

1304964. 

3625. 

.4800-03 

.1169-01 

311927. 

11125. 

.  16C0-03 

.2827-01 

CJi 

CO 


A  EDC-TDR-64-50 


NITROGEN  0 

INPUT 

RUN 

NUMBER  7 

RADIUS 

.50 

1964  NITROGEN  0  DOT 

PROGRAM 

• 

REFERENCE 

DENSITY  0 

.07806410 

LBS/CUBIC  FOOT 

CONDITIONS  BEHIND 

NORMAL 

SHOCK 

T  SEC 

P-2  PSI 

RHO-2  ATM 

T— 2  DEG 

K 

H-2  SOFT/SQSEC 

U-2  FT/SEC 

M-2 

.0000 

12.86374 

.133371 

1768. 

.242935+02 

1000.90. 

.3763 

.0000 

9.19886 

.088630 

1969. 

.274445+02 

1047.39 

.3706 

.0000 

5.53072 

.043846 

.2355. 

.340621+02 

1136.11 

.3639 

.0000 

1.85115 

.009899 

3503. 

.527459+02 

1339.32 

.3243 

INCHES 


.  DATA  PROGRAM 


006A 


CJi 


AEDC-TDR-64-50 


$ 


i 


NITROGEN  Q  INPUT  RUN  NUMBER  3 

1964  NITROGEN  Q  GOT  PROGRAM 

REFERENCE  DENSITY  0.07806410  LBS/CUBIC  FOOT 
RESERVOIR'  CONDITIONS 


RADIUS 


.50 


INCHES 


DATA  PROGRAM 


STAGNATION  CONDITIONS  BEHIND  NORMAL 


T  SEC 

PO  PSI 

RHO-0  ATM 

TO  DEG  K 

HO  SQF  T/SC  SEC 

SOR 

Q-DOT 

PO+PSI 

RHO-0+ ATM 

TO+DEG  K 

SO+/R 

.0000 

16000. 

122.36 

2001 . 

.26225+0 <5 

23.274 

180.00 

14.000 

. 1238183 

2123. 

30.58 

.0000 

16000. 

109.19 

2289,. 

.30262+08 

23.859 

180.00 

10.000 

.0777099 

2390. 

31.48 

.0000 

16000. 

93.20 

2756. 

.36813+08 

24.671 

180,00 

6.000 

. 0390684 

2845. 

32.75 

C71 

Cn 


C06A 


SHOCK 


A  EDC-TDR-64-50 


NITROGEN  Q 

INPUT 

RUN 

NUMBER  8 

RADIUS 

.50 

INCHES 

DATA  PROGRAM  006A 

1964  NITROGEN  Q  DOT 

PROGRAM 

REFERENCE 

DENSITY  0 

.07806410 

L8S/CUBIC  FOOT 

FREE 

STREAM  CCNDITICNS 

MODE  t 

LENGTH  12.000  INCHES 

T  SEC 

P-1  PSI 

RHO-1  ATM 

T— 1  DEG 

K  U— 1  FT/SEC  M-l 

Q-l  PSI 

RE/INCH 

RE/MODEL 

A/ A* 

PQ+/PO  V-l 

•  OOOO 

.0640644 

.0175946 

67.68 

7138.  12.97 

7.545 

261273. 

3135282. 

1951. 

.8750-03  .6828-02 

.0000 

.0411212 

.0108446 

7  C .  4  8 

7678.  13.67 

5.382 

166351. 

1996213. 

2755. 

.6250-03  .9020-02 

.0000 

.0209816 

.0053183 

73.33 

8485.  14.81 

3.223 

86647. 

1039760. 

4658. 

.3750-03  .1354-01 

Ol 


A  EDC-TDR-64-50 


NITROGEN  0  INPUT  RUN  NUMBER  8 

1964  NITROGEN  Q  DOT  PROGRAM 

REFERENCE  DENSITY  0.0780641C  LBS/CUBIC  FOOT 
CONDITIONS  BEHIND  NORMAL  SHOCK 


T  SEC 

P-2  PSI 

RH0-2  ATM 

T,-2  DEG  K 

H-2  SOFT/ SQSEC 

.0000 

12.88849 

.116623 

2086. 

.294030+02 

.0000 

9.21844 

.073237 

2349. 

.339550+02 

.0000 

5.54097 

.036848 

2794. 

.413471+02 

1!  '  -■ 
fe 


RADIUS 


INCHES 


DATA  PROGRAM  D06A 


U-2  FT/SEC 


1076.83 


1136.96 


1224.61 


A  E  DC-T  DR-64-50 


